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DIRECT  MEASUREMENTS  OF  TURBULENT  DISSIPATION 
RATE  USING  DUAL  LASER  DOPPLER  VELOCIMETERS 


INTRODUCTION 

The  leseafch  program  had  two  major  thrusts,  both  of  which  have  been  con^eted.  The  first  was  to 
perfonn  a  detailed  analysis  of  two-point  and  single  point  velocity  measumnents,  obtained  using 
laser  Dealer  velocimetry  (LDV),  in  an  effort  to  give  accurate  turbulent  length  scales  and  die 
turbulent  dissipation  rale.  Four  puUications  reporting  our  findings  regarding  the  calculation  of 
turbulence  sca^  form  LDV  measurements  have  result^  in  the  past 

Benedict,  L.  H.  and  Gould,  R.  D.  (1993)  “Calculation  cf  One-Dimensional  Energy 
Spectra  from  Spatial  Correlation  Measurements  using  LDV,“  Accepted  ASME 1993 
Forum  on  Turbulent  Flows,  Washington,  DC,  June  20-24. 

Gould,  R.  D.,  Nejad,  A.  S.  and  Ahmed,  S.  A.,  (1992a),  “TUrbulent  Length  Scale 
Measurements  in  an  Axisymmetric  Sudden  Expansion  Using  LDV,  The  Fourth 
International  Conference  for  Fluid  Mechanics,  Alexandria,  Egypt,  April  28-30. 

Gould,  R.  D.,  Benedict,  L.  H.,  Nejad,  AS.,  and  Akmcd,  SA.,  (1992b),  “TWo-Point 
Velocity  Correlation  Measurements  in  an  Axisymmetric  Sudden  Expansion  Using 
LDV,"  Proc  of  6th  Ind  Sym  on  Appl  of  Laser  Anemometry  to  Fluid  Mech.,  Lisbon, 

Port. 

Goidd,  R.  D.,  Benedict,  L.  H.,  (1992c)  “A  Comparison  of  Spatial  Correlation  and 
Autocorrelation  Measurements  in  an  Axisymmetric  Sudden  Expansion  Using  LDV," 

13th  Symposium  on  Turbulence,  Uruv.  ^  RoUa-Missouri,  RoUa,  Mo.,  September 
21-23. 


A  summary  of  the  in^rtant  findings  from  this  analysis  is  given  below.  The  second  task  invdved 
fabricating  a  dual  LDV  cqiable  of  making  two-point  and  single-^int  velocity  correlation 
measurements  at  the  ^jplied  Energy  Research  Laboratory  (AERL)  at  North  Carolina  State 
University  (NCSU).  Details  of  this  system  are  also  givoi  below. 

EXPERIMENTAL  APPARATUS 

Two-point  velocity  correlation  measurements  using  two  single  con^ponent  LDV’s  and  single-point 
velocity  correlation  measurements  using  one  single  con^xxient  LDV  were  made  in  an  axisymmetric 
sudden  expansion  air  flow^.  The  geometiy  was  produced  by  joining  a  101.6-mm  (4-in.)  inside 
diameter  entry  pqie  to  a  lS2.4-mm  (6-in)  inside  diameter  clear  acrylic  test  section.  TIk  entry  pipe 
was  3.S-m  long  so  that  a  friUy  develr^red  p^  flow  velocity  profile  existed  at  the  entrance  to  dw 
sudden  expansion.  The  st^  height  for  this  geometiy  was  25.4-mm  (l-ia).  Flat  quartz  windows 
50-mm  X  152-mm  x  3.2-mm  (2  x  6  x  0.125-in)  were  mounted  in  flanges  on  bodi  sides  of  the 
152.4-mm-diameter  test  section  such  that  the  initer  flat  surfaces  were  flush  with  the  inside  diameter 
of  the  test  section. 
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Flow  ^ocities  for  tibe  q;>atial  condadons  were  measured  using  two  TSI  single  confKxient 
dual-beam  LDV  systems,  bodi  operating  in  backscatter  nnode.  Bodi  systems  were  oriented  to 
measure  die  axial  velocity  component  on  the  diaroeter  cf  the  test  section  as  shown  in  Figure  1.  One 
LDV  i^slem  was  adjusted  so  that  the  probe  volume  was  located  at  the  required  axial,  x,  and  radial, 
r,  measurement  locittkm.  Once  this  positkm  was  found,  the  LDV  system  was  locked  m  jdaoe.  The 
514.5-pm  laser  line  from  a  Model  ^2S  Spectra  Riysics  argon  ion  laser  was  used  in  this  system.  A 
Bragg  cell  shifted  the  frequency  of  one  beam  by  40  MHz  causing  die  fringes  to  move  in  the 
downstream  direction.  A  seccmd  LDV  system  (TSI  Model  9277  19(Mnm  fiber  opdc  i»obe), 
mounted  on  a  {»ecision  xyz  positioning  tai^  with  resdution  of  ±  2.5pm  in  eadi  axis,  was  located 
on  die  opposite  side  of  the  test  sectkm  (see  Figure  1).  The  488-pm  laser  line  frtm  a  Model  165 
Spectra  Pt^sics  argon  ion  laser  was  used  in  tfali  system.  A  frequency  diifi  of  40  MHz  was  used 
causing  the  fringes  of  diis  syston  to  move  in  die  iqistreon  diiecdrm.  Bodi  IDV  systons  employed 
3.75x  beam  expansion  op^  and  gave  probe  volumes  qiproximately  60-pm  in  dianeter  aid 
4S0-pm  in  lengfo. 

A  20-pm-diameter  pinhole  mounted  on  a  fixture  supported  on  a  spare  test  section  window  was  used 
to  find  the  position  where  bodi  laser  beam  probe  volumes  initially  overlapped.  This  fixture  was  used 
prior  to  each  test  sequence,  dius  ensuring  dia  both  probe  vdumes  overlaf^ied  a  the  zero  sqiaraion 
distance  point.  Specially  designed  beam  blocks  were  fabricated  to  block  reflections  from  die  LDV 
focusing  lens  (diey  faoc^  one  another)  and  from  test  section  windows.  Narrow  bandpass  filters  were 
placed  in  front  of  each  photomultiplier  tube  to  eliminate  cross-talk  between  the  two  charmek. 

IVo  TSI  Model  1990C  counter  processors  interfaced  to  a  custom  built  coincidence  timing  unit  were 
used  in  the  data  collection  and  processing  system.  Hi^  and  low  pass  filters  were  set  to  10  MHz  and 
50  MHz,  mpectively,  for  the  stat^nary  LDV  system,  and  20  Nfiiz  and  100  MHz,  reflectively,  for 
the  fiber  optic  LDV  system  Both  processors  were  set  to  sanqde  oontinuou^  (midtqde 
measurements  per  burst)  count  16  fringes  and  use  a  1%  conoparator.  A  hardware  coincideiit  window 
was  set  at  20-ps  for  all  of  the  tests.  Data  (two  velocities  and  die  running  tnite  for  eadi  realization) 
were  transferr^  through  two  DMA  ports  to  a  MicroVax  miniconqiuter  later  uploaded  to  a  VAX 
8650  for  analysis. 

The  flow  field  was  seeded  using  titanium  dioxide  Cl^)  particles  generated  by  reacting  dry 
titanium  tetra-dilorkle  (1104)  w^  the  moist  shop  air.  Oaig  et  al.^  measured  the  particle  sizes 
generated  by  this  device  and  found  diat  they  were  fairly  unifonn  and  in  the  0.2  -  1-pm-diameter 
range.  Data  validation  rates  varied  between  5000  and  500  per  second  on  each  counter  processor  and 
deperuied  mainly  on  how  well  die  chemical  reaction  proceeded.  This  seemed  to  be  very  sensitive  to 
sh^  air  temperature  and  relative  humidity.  Coincident  data  validation  rates  tanged  from  1000  to  50 
measumnents  per  second. 

EXPERIMENTAL  PROCEDURE 

An  flow  conditions  were  maintained  at  near  constant  values  throughout  die  testing  procedure.  The 
inlet  centeriine  velocity,  Uc/,  was  maintained  at  18.0-m/s  ±  0.1-in/s  (59-fr/s  ±  0.3-fos)  giving  Rcd  = 
114,000  based  on  coiterline  velocity  and  inlet  diameter.  Spatial  correlation  statistics  and  histograms 
were  formed  by  using  5000  individual  realizaticais  for  each  velocity  channd  at  each  measurement 
point.  Autocorrelations  were  formed  by  using  50,000  individual  realizations  from  die  stationary 
LDV  system  widi  data  validation  rates  between  25,000  and  1000  per  second,  fri  confuting 
statistical  parameters,  a  two  step  process  suggested  by  Meyers^  was  used  to  eliminate  noise  fr^  the 
data.  All  velocity  measurements  deviating  more  than  3  standard  deviations  from  die  mean  are  dius 
eliminated.  For  aprcperly  operating  LDV  system  less  dian  1%  of  the  data  should  be  discarded.  In 
addition,  Srikantahdi  and  Coleman*  have  found  diat  as  long  as  the  total  number  of  removed  points 
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Figure  1.  Dual  LDV*s  in  an  axisymmetric  sudden  expansion 


Figure  2.  Autocorrelation  measurements. 
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is  within  1%  of  the  total  satire  sue,  the  effect  of  the  removal  of  noise  tm  the  autoc(»relation 
function  is  negligiUle. 

EXPERIMENTAL  RESULTS 

Autoconelations  and  Mkaoscaks 

l\vo-point  (spatial)  velocity  correlation  measurements  and  single  point  (auto)  velocity  conelation 
measurements  were  made  at  three  locations  in  the  axisymmetric  sudden  expansion  flowfield.  These 
results,  including  the  calculatkm  of  integral  lengdi  scales  and  Taylor  microscales,  were  presented  by 
Oould  and  Benedict^.  In  this  earlier  piper  mention  was  macle  of  the  difficulty  in  accurately 
determining  die  Taylor  microscale  from  eidier  die  discrete  qntial  or  autocorrelation  frmctioiis.  This 
difficulty  is  addressed  in  more  detail  tteie.  Additiraially,  the  problem  of  calculating 
tme-duTHmsional  energy  spectra  from  spatial  and  autocorrelation  hinctions  is  discussed. 

In  Figure  2,  typical  autocorrelation  functions  measured  at  6  step  heights  downstream  (x/H  =  6)  and 
2  step  heig^  (r/H  =  2)  from  die  axis  are  presented.  The  autocorrdation  can  be  transformed  to  die 
spatial  domain  by  use  of  Taylor’s  hypothesis  (i.e.,  x=Vt )  which  was  shown  to  be  an  accurate 
transfoimation  for  this  flow  fields  Discrete  autocorrelation  measurements  were  made  using  die 
slotting  tedinique  described  by  Jones^  and  Mayo,  et  al.^  The  time  lag  axis  was  divided  into  bins  of 
equal  width  arid  the  exact  lag  products  of  all  points  up  to  a  specffied  maximum  lag  time  were 
accumulated  in  appropriate  bins.  The  average  of  all  the  auto-products  falling  in  eadi  bin  was 
assumed  to  be  the  value  of  the  discrete  autocorrelation  frmcdon,  Le: 


at  the  mit^iirt  of  the  bin.  The  data  was  filtered  first  using  the  previously  described  mediod  to 
eliminate  velocity  measurmiaits  more  dian  3  standard  deviations  from  the  mean.  The  zero-lag 
{uoducts  were  not  included  to  avoid  the  spiking  effect  at  the  origin  (t  =  0)  due  to  squaring  of 
uncotrelated  noise'^ .  The  discrete  autocoirdatims,  shown  in  Figure  2,  were  obtained  using  a  slot 
width  of  SO-ps  with  a  maximum  lag  time  of  3-ms.  Note  that  there  is  still  some  scatter  in  the  data 
even  diough  50,000  sarrples  were  used  to  build  the  autocorrelation  functions.  A  mean  data 
validation  rate  of  qjptoximately  12,000  per  second  was  used  here.  This  scatter  is  statistical  in 
nature  and  is  due  to  too  few  lag  products  in  each  slot  r^ch  causes  some  variance  in  die  data  (each 
bin  contained  ipproximately  30,000  lag  products).  The  error  bars  shown  in  tfiis  figure  uKlicate  the 
variances  of  ^  calculated  mean  values  of  the  autocorrelation  function.  These  statistical 
uncertainties  in  the  mean  td  the  quantity  m'ii'  are  estimated  using  die  sample  size  and  standard 
deviation  of  die  sairple  in  each  slo^ .  Interestingly,  although  zero-lag  products  were  not  included  in 
dte  first  slot  of  the  discrete  autocorrelation  function,  a  severe  spike  results  at  die  origin  (bottom 
curve).  This  spike  is  an  unnatural  characteristic  of  the  autocorrelation  function  and  {nohibits  die 
estimatkm  of  ^  Taylor  micro-scale  (diss^ation  scale). 

One  pqiular,  but  peiluqis  questionable,  way  around  this  problem  is  simply  to  ignore  the  first  slot 
value  md  normalize  the  autocorrelation  with  respect  to  the  second  slot.  Ano^r  qjproach  was 
considered,  however,  based  <»i  the  belief  that  die  spiking  effect  was  caused  by  the  continuous 
sanqiling  of  the  LDV  system.  Ctmtinuous  sanqiling  allows  for  mult^le  measurements  to  be  made 
on  a  single  Dop[der  bt^.  In  such  a  situation,  a  slow  moving  particle  can  be  sanqiled  more  than 
once  as  it  passes  through  the  probe  vtdume.  Such  measumnents  are  inherendy  hi^y  correlated 
wrdi  one  another  as  ttey  come  from  the  same  particle  at  close  proximity  in  time.  A  simple 


4 


coiqputer  algorithm  was  en:q>loyed  to  remove  imdtqx  velocity  measurements  resulting  from 
passage  of  single  scattering  source  dirough  the  probe  vdume  from  the  data  set.  Basically,  the  data 
set  was  search^  for  successive  measurements  vriiich  had  velocity-time  between  data  im)^cts  less 
than  the  probe  volume  diameter.  If  this  occurred,  these  measurements  were  assumed  to  have  come 
from  die  same  particle.  After  eliminating  these  redundant  measurements,  a  new  autoc<»relatkxi 
function  was  calculated  and  is  also  shown  in  Figure  2  (t<^  curve).  One  may  observe  that  this  new 
autocorrelatitm  fiinctimi  dis|days  the  expected  paraboUc  behavior  near  the  origin  without  any 
spiking. 

From  such  an  autocorrelation  function  the  Taylor  microscale  may  be  easily  estimated  by  fitting  an 
osculating  parabola  to  the  spatial  correlation  function  betweoi  ^  first  measured  two  points.  The 
discrete  spatial  correlation  ftinction  is  obtained  by  transforming  the  discrete  autocorrelation  function 
via  Taylor’s  hypothesis  (le.,  ’Kf^V  XeY  Thus  the  longitudinal  microscale  is  obtained  from: 

where  X  rejxesents  die  separatirm  distance  between  the  first  two  measurements.  For  the 
autocorrelation  function  shown  as  the  top  curve  in  Figure  2,  where  V  -  8.5  m/s,  Xy  was  found  to 
equal  3.3-mm.  In  previous  wotk^  diis  estimate  had  been  considered  to  be  on  the  order  of  6-mm, 
however,  an  earlier  experimental  investigation^  estimated  the  dissqiation  at  this  location  in  the  flow 
field  to  be  t^qiroximately  1340-tnVs^  by  using  turbulent  kinetic  energy  balance.  Using  an  isotropic 
turbulence  assuirption,  the  microscale  can  be  estimated  froiiu 
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The  microscale  estimate  based  rei  a  turbulent  kinetic  energy  balance  and  the  assumption  of  isotropy 
is  then  calculated  as  approximately  2.3-mm.  As  the  dissipation  rates  obtained  from  the  turbulent 
kinetic  energy  balance  are  believed  to  be  more  accurate  than  die  first  microscale  estimates  fiom  the 
autocorrelation  functions  presented  by  Gould  and  Benedict^  the  newer  estimate  of  the  microscale 
from  these  functions  is  considered  to  be  an  inprovement  ovn*  the  old  one. 


As  will  be  shown  later,  two-point  patial  correlations  seem  to  suffer  from  probe  volume  integration 
effects  such  that  near  dw  origin  (Ax  =  0)  no  parabolic  region  can  be  found  and  any  estimation  of  the 
microscale  must  be  considered  a  crude  q>ptoximat^of^  Until  these  probe  volume  integration  effects 
can  be  reduced,  the  transformed  autocorrelation  function  (using  Taylor’s  hypothesis)  riiould  be 
considered  as  the  mediod  of  choice  fra:  estimating  the  micro-scale  from  LDV  measurements. 


Some  guidelines  are  thus  in  order  as  to  how  to  best  pursue  this  goal.  First,  one  should  always  be 
aware  that  Taylor’s  hypothesis  is  limited  by  turbulence  inteisity.  It  should  not  be  pplied  in  flows 
with  turbulence  intensities  greater  than  about  20%.  Work  by  Gould  and  Benedict^  and  Cettedese’^*’ 
suggests,  however,  that  oftentimes  Taylor’s  hypothesis  is  t^proximately  valid  for  much  higher 
turbulence  intensities.  Another  point  to  make  is  that  a  very  high  mean  data  validation  rate  is  often 
necessary  to  resolve  die  small  scales.  Previous  work  such  as  that  of  Cenedese^^,  Morton  and 
Clark^^  Absil'^,  and  Fraser  et  al.^  make  use  of  relatively  low  Reynolds  number  flows  to  simplify 
the  measuring  process.  In  the  current  work,  a  Reyncdds  nuniber  based  on  inlet  diameter  of 
ajproximately  114,000  was  adiieved  in  the  axisymmetric  sudden  expansion  flow.  For  this  flow 
condition  additional  studies,  not  repotted  here,  suggest  that  pproximately  100,000  sanples  need  to 
be  recorded  at  data  rates  ppto«d^  20,000  Hz  in  order  to  accurately  resolve  the  small  scales. 
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Oftentimes,  however,  such  high  data  rates  are  not  achievable.  Additkmally,  single  measurement  per 
burst  mode  is  preferred  over  continuous  sani{ding  in  (»der  to  prevent  biasing  of  die  autoconelation 
function  as  de^bed  above.  As  was  deme  in  this  work,  zero-lag  products  should  not  be  kiduded  in 
the  first  slot  of  die  discrete  autoonielation  function. 

Spatial  Correlations  and  1-D  Energy  Spectra 

As  mentioned  above,  ^atial  conelations  suffer  from  a  probe  volume  integration  effn^t  at  amall 
separation  distance.  This  effect  tends  to  obscure  correlation  information  pertaining  to  the  microscale 
as  even  whm  die  two  probe  volumes  overliq)  (i.e..  Ax  =  0),  their  lengdis  allow  frv  a  finite  sqiaration 
distfflice^  If  the  probe  volume  length  is  of  die  same  size  as  the  par^lic  region  of  the  c<Hielatk»i 
function,  then  this  parabolic  section  can  quite  possibly  be  obsoired.  As  flow  Reynolds  rnimbo’ 
increases,  such  a  situation  is  more  likely  to  hsi^ien,  aikl  thus  die  determination  of  a  microscale 
dirough  a  two-point  ^atial  correlation  function  becomes  an  iq^iroximate  analysis. 

Figure  3  shows  the  longitudinal  spatial  correlatkm  coefficients  at  x/H  =  10  and  r/H  =  2  in  the 
axisymmetric  sudden  expansion  flow  obtained  usmg  two-point  correlation  measuronents  and  also 
shows  an  exponential  fit  to  the  data.  The  lack  of  a  parabolic  region  in  the  data  near  zero  ^latial 
separation  should  be  pointed  out.  Also,  it  can  be  seen  that  die  ^atial  correlation  function  do^  not 
reach  unity  at  zero  spatial  separation.  Bodi  of  these  diaracteristics  may  be  attributed  to  die  probe 
volume  integratiem  effect  mentioned  above.  Obviously,  diere  is  litde  basis  for  fitting  a  parabola  to 
such  a  correlatirm  function  near  zero  separation,  whidi,  in  turn,  increases  the  difficulty  in  making  an 
estimation  of  the  microscale. 

It  is  also  possible  to  make  an  estimation  of  the  microscale  from  the  one-dimensional  energy 
spectrum  (fefined  as  the  Fourier  transform  of  the  spatial  catrelaticai  function: 

Ey{k)  =  R^^{x^)cosikx^) 

where  TT^  is  the  variance  of  die  velocity,  77  is  the  mean  velocity,  x  is  the  separation  distance  and  k  is 
the  wavenumber.  The  subscrqit  "1"  denotes  that  quantities  are  measured  in  die  axial  direction  only. 
The  microscale  can  be  obtaiiKd  by  integration  of  the  (xie-dimensicmal  energy  spectrum 


where  A/  denotes  die  integral  length  scale  in  die  axial  direction. 

From  die  correlation  function  shown  in  Figure  3,  the  1-D  energy  spectrum  may  be  arrived  at  by  a 
discrete  numerical  integration  process  involving  the  cosine  transform  of  the  correlation  function  as 
stated  mathematically  above.  This  process  may  be  carried  out  on  the  discrete  correlation  data  points 
themselves  if  care  is  exercised.  If  the  discrete  cosine  transform  is  tpplied  to  the  discrete  correlation 
data  shown  in  Figure  3,  an  unrealistic  energy  spectrum  as  shown  in  Figure  4  results.  The  solid  line 
in  Figure  4  shows  the  dieoretical  energy  spectrum  corresponding  to  the  exponential  fit  to  the  spatial 
cturelaticm  data. 
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Figure  3.  Longitudinal  spatial  correlation  measurements. 
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Figure  4.  I-D  energy  spectrum. 
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Figure  5.  1-D  energy  spectrum. 


Wavenumber,  k  (lAn) 


Figure  6.  Lateral  spatial  correlation  measurements. 
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This  misrepresentation  of  die  1-D  eneigy  spectrum  can  be  traced  to  two  problem:  (1)  a  non-zero 
value  for  the  correlation  coefficient  at  maximum  separation  distance  (i.e.,  the  correlation  function 
does  not  approach  zero  closely  enough)  and,  (2)  insufficient  spatial  resolution  of  the  cosine  function 
within  the  transform  at  high  wavenumbers  (i.e.,  discrete  con^ation  coefficients  not  spaced  closely 
enough).  The  first  problem  is  usually  solved  by  inlying  a  window  function  to  the  correlation 
coefficient  such  that  the  correlation  coefficient  is  forc^  to  zero  when  Ax  becomes  larger  dian  a 
prescribed  value. 

To  study  the  problem  of  nonzero  correlation  coefficient  at  maximum  separation  distance,  die 
exprmential  fit  to  die  data  was  integrated  analytically  using  various  maximum  sqparation  distances. 
Analytic  integration  allowed  the  truncation  |»oblem  to  be  isolated  fix»n  die  spatial  resolution 
inaccuracies.  Figure  5  shows  the  effect  produced  by  integration  over  a  distant^  equivalent  to  the 
maximum  s^aradon  present  in  the  data  set  shown  in  Figure  3  (i.e.,  Axmox  =  88-mm).  Since  the 
correlation  fimction  did  not  reach  zero  at  this  maximum  separation  distance  (the  LO  Vs  could  not  be 
separated  furdier),  a  tnmcation  effect  in  the  Fourier  transform  produced  the  oscillatory  behavior  in 
the  oiergy  spectrum  as  seen  in  Figure  S.  This  proUem  can  be  eliminared  through  the  use  of  a 
smoothing  window  which  effectively  drives  the  correlation  coefficient  to  zero  at  maximum 
separation. 

To  study  the  resolution  problem,  the  correlation  function  was  modeled  widi  an  exponential  fit  with 
various  discrete  separation  distances  in  the  data.  The  Hanning  window  was  used  to  eliminate  the 
problem  of  truncation  of  the  correlation  curve.  Figure  6  shows  die  discrete  cosine  transform  results 
for  an  exprmential  correlation  coefficioit  when  35  evenly  spaced  points,  out  to  a  separation  distance 
of  88-mm,  were  used.  Notice  the  large  spike  in  the  one-dim^ional  energy  spectrum  at  high 
wavenumber.  This  corresponds  to  the  wavenumber  at  which  the  cosine  function  in  the  Fourier 
transform  is  no  longer  prqieiiy  resolved.  Proper  resolution  of  the  cosine  function  requires  the 
spacing  of  the  correlation  data.  Ax,  to  be  less  than  itiTk. 

The  above  exanqiles  indicate  that  fairly  siirqile  operatkms,  such  as  die  Fourier  transform  of  a 
correlation  function,  can  give  highly  erratic  resets  if  used  improperly.  In  order  to  properiy  rqiresent 
the  one-dimensional  eneigy  spectrum  fiom  spatial  corrdation  fiinaions,  the  correlation  functitm 
must  be  driven  to  zero  by  a  window  function  if  it  does  not  go  to  zero  naturally  at  the  maximum 
separation  distance.  Note  diat  there  will  always  be  apiactical  maximum  sqiaration  distance  limiied 
by  the  experimental  equqnnent.  Moreover,  dre  spacing  of  the  discrete  data  points  must  be  less  than 
rb2k,  where  k  is  die  maximum  wavenumber  of  interest. 

The  maximum  wavenumber  for  which  the  energy  spectrum  may  be  calculated  can  be  increased  by 
curve  fitting  experimental  data  and  interpolating  intermediate  v^ues  of  the  correlation  coefficient  - 
while  remaining  conscious  of  the  n/2k  guideline,  ft  is  well  known  that  the  exponential  fit  shown  in 
Figure  3  has  a  Fourier  cosine  transform  which  may  be  integrated  analytically^^.  This  transform  is 
represented  as  either  a  solid  or  dashed  curve  in  Figures  4-6.  ft  should  be  not^  diat  the  exponential 
representation  of  die  spatial  correlation  function  produces  a  1-D  energy  spectrum  with  a  slope  of  -2 
on  a  log-log  axis  plot.  Such  a  function  has  a  non-finite  integral.  The  Taylor  microscale,  therefore, 
caimot  be  estimated  by  integration  of  the  1-D  eneigy  spectrum  whidi  is  calculated  using  a  pure 
e^xMiential  representation  of  the  spatial  correlaticm  fimction. 

In  order  to  more  properly  represent  the  1-D  eneigy  spectrum  a  multi-piece  fit  to  either  the  ^atial  or 
autocorrelation  function  is  necessary.  Since  the  hi|^  fiequency  portion  of  the  spectrum  is  dominated 
by  bdiavior  near  zero  separation  in  the  correlatitm  fimctions  and  typical  correlation  functions 
possess  a  parabolic  region  near  zero  separation  (not  characterized  by  die  exponential  fit  shown  in 
Figure  3),  ap  attend  was  made  to  mate  a  parabolic  curve  near  zero  separation  to  an  exponential 
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curve  (whidi  represents  die  measurements  well  at  large  separation  distances).  Some  results  qipear 
in  Figure  7,  wldch  ^ow  die  effects  on  die  1-D  energy  ^lectrum  of  mating  pandiolic  regions 
iqnesenting  differing  microscales  to  die  exponential  curve.  As  eiqiected,  the  parabdic  region  is 
resptmsible  for  driving  the  energy  spectrum  to  zoo  at  high  fiteqiiency  and  that  die  huger  die 
microscale,  the  lower  ^  wavenumber  where  this  luqipens.  This  inakes  Mnse  from  the  stand^int 
that  a  larger  microscale  in^lies  more  correlation  which,  in  turn,  implies  less  turbulence  intensity  (in 
the  limit  laminar  flows  are  perfeedy  correlated)  and  hence  less  energy.  Interestingly,  die  steep 
descent  of  the  1-D  energy  spectrum  at  high  frequency  mimics  die  (-7)  8l<^  predicted  by 
Kolmogoroff  for  homogeneous  turbulence. 

These  energy  spectra,  which  result  when  representing  die  qiatial  conelation  codficioit  with 
two-piece  fits,  may  indeed  be  integrated  using  Eq.  (5).  The  microscales  returned  die  inregration 
process,  however,  are  30  to  40  percent  higher  dian  the  original  values  assumed  to  determine  the 
exact  stupe  of  ^  parabolic  fit.  Obviously,  an  actual  correlation  function  cannot  usually  be 
represent^  exaedy  by  mating  a  parabolic  and  an  exponential  curve;  so  one  would  not  expect  sudi  a 
fit  to  return  die  exact  value  of  the  microscale  when  determined  by  numerical  integration  of  the  1-D 
tnagy  spectrum.  Microscales  calculated  by  numerical  integration  of  the  1-D  energy  spectrum 
would  be  reasonably  close,  however,  to  the  actual  values  and  would  provide  valuable  insight  into 
turbulent  structure. 

Figure  9  is  a  similar  plot  of  die  1-D  energy  spectrum,  fri  this  case,  the  1-D  energy  spectrum  using 
the  paraboUc  fit  near  zero-separadm  is  compared  to  that  using  a  third  order  polynomial  fit.  In  both 
cases,  these  curves  are  mated  to  the  same  exponential  curve.  Note  diat  die  Aird  Mder  fit  drives  die 
1-D  energy  spectrum  to  zero  at  a  lower  wavenundier  when  conpared  to  a  parabolic  fit.  Both  fits 
were  derived  from  an  assumed  miaoscale  of  lO-tnm.  Numerical  integration  of  the  1-D  energy 
spectnim  (using  Eq.  S)  returned  a  miaoscale  value  of  12.6-mm  for  the  diird  order  fit  as  opposed  to 
14.2-nim  for  the  parabolic  fit.  Once  again  this  trend  makes  sense  as  the  third  order  fit  rolls  off  more 
slowly  than  a  parabolic  fit  giving  higher  coreladon  values  near  zero  sqiaration  and  ther^  less 
turbuloit  intensity  and  less  energy.  Evidendy  a  fourdi  or  fifdi  order  fit  would  return  through 
numerical  integration  of  die  1-D  enogy  spectrum  almost  die  exaa  value  of  the  microscale  diat  h^ 
been  used  to  create  the  fit.  It  is  recommmded  then  that  when  using  LDV,  microscales  be  determined 
by  the  following  procedure:  (1)  measure  the  autocorrelation  function,  (2)  fit  this  function  with  a 
higher  order  polynomial  near  zero-separation  and  an  erponential  fitnetion  at  larga  sqiaration 
distances,  (3)  calculate  the  1-D  enagy  spectrum  from  the  curve  fit  to  the  conelation  function,  md 
(4)  numerically  integrate  the  resulting  1-D  energy  pectnim  to  determine  die  microscale.  When 
making  the  curve  &  to  the  experimental  data  and  interpolating  between  data  points,  pacing 
requirements  and  window  function  requirements  should  be  observed  as  moitioned  above. 

One  last  point  to  mention  concerning  the  curve  fitting  of  correlation  functions  is  that  lidien  merging 
two  different  curve  fits,  continuity  of  slope  is  of  utmost  inportance.  Any  discontinuity  in  slope  will 
produce  ringing  at  hi^  wavenumber  in  die  calculated  1-D  energy  pectrom.  This  is  because  the 
Fourier  cosine  transform  of  the  correlation  fiuiction  is  the  integral  of  the  correlatioi  function 
multplied  a  cosine  funaion  which  in  efifea  is  a  weakly  danped  cosine  function.  As 
wavenumber  increases  such  a  function  has  an  integral  which  rpidly  pproaches  zero,  which  is 
profoundly  obvious  when  the  energy  spectrum  is  plotted  on  normal  (not  log)  axes.  A  small 
discontinuity  in  the  sl<pe  of  die  correlation  functioi  thus  changes  die  way  the  integral  pproaches 
zero.  But  for  each  new  wavenumber  this  hppens  in  a  different  way  since  die  discontinuity  occurs 
eadi  time  at  a  different  location  on  die  damped  cosine  function.  Ringing  in  die  1-D  energy 
spectrum  is  a  direa  and  inevitable  coisequence  of  slope  discontinuity. 
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Figure  7.  1-0  energy  spectrum(fits  to  correlation  function). 


Figure  8.  1-D  energy  spectrum(fits  to  correlation  function). 
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Figure  10.  Single-component  fiber  optically  coupled  LDV  system. 
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DEVELOPMENT  OF  DUAL  LDV  AT  NCSU 

Optical  compmieiits  including,  fiber  optic  input  and  ou^t  couplers,  polarization  preserving  optical 
fiber,  banc^^^  laser  line  filters,  and  an  acoustic  optic  modubttor  (Bragg  cell)  were  purcha^  to 
modify  an  existing  LDV  system  so  diat  two-point  wlocity  conelatkHi  measurements  can  be  made  at 
NCSU.  A  sdwmatic  diagram  of  the  complete  dual  LDV  system  is  shown  in  Figures  9  tnd  10.  The 
LDV  system  shown  in  Figure  9  is  a  two-color,  two-conqtonent.  It  has  been  modified,  by  adding 
beam  splitting  prisms  and  fiber  optic  ii^xit  coulters,  to  provide  laser  light  for  the  second  LDV 
system  which  is  required  to  make  two-point  velocity  measurements.  71113  seamd  system  is  a 
single-col<v  single-component  LDV  system.  Note  that  two  fiber  optic  input  coujders  have  been 
adM  to  the  two-color  system.  This  allows  the  flexibility  to  choose  either  color  (blue  or  green)  fx 
the  single  component  system.  Both  systems  are  mounted  on  three-axis  tables  so  that  the  probe 
volumes  can  be  placed  at  user  selected  locations.  Presently,  both  systems  operate  in  forwad  scatter 
collection  mode;  however,  modificaions  are  being  made  to  the  systems  so  that  either  forward 
scatter  or  90  degree  coUectitm  angle  can  be  used.  In  addition  to  die  optical  conqionents,  the  LDV 
data  acquisition  interface  has  been  upgraded  so  diat  up  to  three  channeb  of  velocity  data  can  be 
acquired  simultaneously.  It  is  envisioned  that  diis  {pparatus  will  be  used  to  make  many 
stae-of-the-ait  turbulence  measurements  in  the  future.  Nfeasurements  in  combustion  environments 
are  of  particular  interest 

SUMMARY 

It  has  been  shown  that  spiking  in  autocorrelation  measurements  obtained  using  LDV  is  a  prevalent 
phenomena  that  may  be  caused  by  biases  induced  by  continuous  samjding.  This  qnking  should  be 
distingui^ied  fiom  that  attributed  to  inclusion  of  zero-lag  terms  in  the  slotting  tedinique.  The 
piking  caused  by  a  continuous  san^ling  bias  may  be  overcome  by  filtering  die  data  to  remove  the 
mult^le  sampled  velocity  realizations.  Spatial  cxrelation  measurements  have  been  shown  on  (he 
other  hand  to  sulGfer  fiom  probe  volume  integration  effects  which  cannot  be  remedied  once  the 
measurements  have  been  made.  These  probe  volume  integration  effects  may  obscure  the  dassic 
parabolic  region  near  zero-separation  in  spatial  ctarelation  functknis,  which  degrade  micro-scale 
estimations. 

Calculation  of  1-D  eneigy  spectra  fiom  cozrelatkm  functions  for  purposes  of  calculating  the  Taylor 
mkroscide  is  a  viaUe  alternative  to  fitting  a  parabda  to  either  ^atid  or  autocondation  functim. 
It  has  been  shown,  however,  diat  even  a  faiily  sinqde  operation  such  as  the  Fourier  transform  of  a 
correlation  function  can  give  highly  erratic  results  if  used  inqnopetly.  In  order  to  properly  represent 
eneigy  spectra  fiom  correlation  measurements  or  curve  fits  to  those  measurements,  the  correlation 
function  mur'  be  driven  to  zero  by  a  window  function  and  the  spaciiig  of  die  discrete  data  points 
must  be  less  than  n/2k,  where  A:  is  the  maximum  wavenumber  of  interest.  For  piece-wise  curve  fits, 
continuity  of  slope  must  be  maintained  when  merging  die  different  pieces  of  a  given  curve  fit. 
Lastly,  an  in-house  dual  L  HV  system  cipable  of  making  getreral  two-point  velocity  measuremenis 
has  been  fabricated  and  will  be  used  to  further  our  oncferstandiiig  of  turbulence  and  die  cond^jstion 
process. 
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